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Introduction 

In this paper, we consider a multidimensional mixed stochastic differential equation 
of the form 

X{t) =X + £ (a(s,X(s))ds + b(s,X(s))dW{s)+c(s,X(s))dZ(s) S j, (1) 

A. Melnikov 

Department of Mathematical and Statistical Sciences, University of Alberta, 632 Central Academic Build- 
ing, Edmonton, AB T6G 2G1, Canada 
E-mail: melnikov @ualberta.ca 

Yu. Mishura 

Kyiv National Taras Shevchenko University, Faculty of Mechanics and Mathematics, Department of Prob- 
ability, Statistics and Actuarial Mathematics, Volodymyrska 64, 01601 Kyiv, Ukraine 
Tel.: +380-44-259-03-92 
Fax: +380-44-259-03-92 
E-mail: myus@univ.kiev.ua 

G. Shevchenko 

Kyiv National Taras Shevchenko University, Faculty of Mechanics and Mathematics, Department of Prob- 
ability, Statistics and Actuarial Mathematics, Volodymyrska 64, 01601 Kyiv, Ukraine 
Tel.: +380-44-259-03-92 
Fax: +380-44-259-03-92 
E-mail: zhora@univ.kiev.ua 



2 



Alexander Melnikov et al. 



where W is a standard Wiener process, and Z is an adapted process, which is almost 
surely Holder continuous with exponent y > 1/2. 

The strongest motivation to study such mixed equations comes from financial 
modeling. The observations of stock prices processes suggest that they are not self- 
similar: on a larger time scale (months or years) these processes are smoother and 
have a longer memory than on a smaller time scale (hours or days). One reason for 
this is that the random noise in the market is a sum of a more irregular "trading" noise 
rendering irrationality of the stock exchange, and a more regular "fundamental" noise 
rendering a current economical situation. The first random noise component prevails, 
especially for illiquid instruments, in a shorter time periods (days and hours). The 
second one takes some time to propagate and becomes essential in a long run, clearly 
exhibiting a long memory. Such phenomena can be modeled by a sum of a Wiener 
process W and a fractional Brownian motion B H with the Hurst parameter H > 1/2. 
The behavior of this process on a smaller scale is mainly influenced by independent 
increments of W and its irregularity, while on the larger scale the long memory of 
B H dictates the evolution of the process. As a result, the mixed model describes the 
stock price behavior in a better way. Let us note that the long memory effect in the 
financial markets and application of the models involving fractional Brownian motion 
and mixed Brownian-fractional Brownian motion were studied in many papers (see 
E), 0, 0, El, (23) and references therein). 

The existence and uniqueness of a solution to ([TJ was established in ifTTI . ifTTl . 
|[T9l and ll20l under different assumptions, and the most general results are obtained 
in IfTTl and [20 1 . Besides this, in ll20l a limit theorem for mixed stochastic differential 
equations was established, and we apply this result here. 

In this paper we study a stochastic viability property of the solution to equation 
([TJ and the applications of this property. A stochastic process X is called viable in a 
non-random set D, if starting at Xq € D, the process stays in D almost surely. For Ito 
stochastic differential equations such property was studied, for example, in UJ, 0, 
0, ED, HSl, (25]. Recently (see flU), a viability result was proved for stochastic 
differential equations with fractional Brownian motion. Studying a viability prop- 
erty of solution of equation (QJ, we deduce the path-wise comparison theorems as an 
important application. The first comparison theorems for Ito stochastic differential 
equations were obtained in [12|, [22 1, [24]. In papers lfl6l and ITSl multidimensional 
pathwise comparison theorems were obtained too (in the first one, for processes with 
jumps). Here we prove a viability result for ([JJ) and deduce the results on positivity 
and pathwise comparison for solutions of such equations. 

The paper is organized as follows. Section 1 contains some preliminaries on gen- 
eralized pathwise stochastic integration. It contains also some results about existence 
and uniqueness of solution to ([TJ as well as a limit theorem for for solutions of such 
equations. The main results on viability, positivity and pathwise comparison of so- 
lutions to mixed stochastic differential equations are presented in Section 2. Section 
3 gives some applications of the results of Section 2 to the comparison of option 
prices in the market models with long-range dependence. More precisely, we consider 
here pure fractional and mixed Brownian-fractional-Brownian Cox-Ross-Ingersoll 
(mBfBm-CIR) models. The existence and uniqueness result for the corresponding 
stochastic differential equations is proved. The comparison theorem is applied to 
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obtain an upper bound for the price of a European-type option on the interest rate 
presented by the mBfBm-CIR model. We also give a viability and comparison re- 
sult for Ito stochastic differential equations with random coefficients. The proofs of 
these results are very similar to those for equations with non-random coefficients, so 
we put them to Appendix, where we also discuss some multidimensional pathwise 
comparison results. 

1 Preliminaries 

Let (X2, J?,F, P) be a complete probability space with a filtration F = {J?f,f > 0} 
satisfying usual assumptions. Let W — { \W\ (?),-•• , W m (t)),t > 0} be a standard m- 
dimensional F-Wiener process andZ= { [Z\ (?), . . . ,Z r (f)),f > 0} be an r-dimensional 
F-adapted /-Holder continuous process on this probability space. 

We consider the following mixed stochastic differential equation in W 1 : 

rt , m r 

X{t)=X + / (a{s,X(s))ds+Yb k {s,X(s))dW k (s) + Ycj{s,X(s))dZj{s)), 

te [0,r], 

(2) 

where coefficients a,^,c,-: [0,T] x M. d — > M. d are jointly continuous. For brevity, we 
will use notation 

m 

b{s,X{s))dW{s) := £ h(s,X(s))dW k (s) 

k=l 

and 

c(s,X(s))dZ(s) :=£cj(s,X(s))dZj(s). 

We note that in equation d2l, the integral w.r.t. Wiener process is defined as the 
standard Ito integral, and the integral w.r.t. fBm is pathwise generalized Lebesgue- 
Stieltjes integral, whose definition is given below. 

1.1 Generalized Lebesgue-Stieltjes integral 

Consider two continuous functions f,g £ C(M. + ). For a £ (0, 1) and < a < b define 
fractional derivatives 

Note that the latter notation is slightly different from the one given in [21]. This 
simplifies the notation below. If f,g are such that D" + f £ L p [a,b], D l b Z a g € Lq[a,b] 
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for some p E (1,1 /a), q = p/(p— 1), we can define the generalized (fractional) 
Lebesgue-Stieltjes integral as 



f(x)dg{x) = e'™ / {D a a+ f) (x) (D)- a g) (x)dx. 



(3) 



As Z is almost surely /-Holder continuous, it is easy to see that for any a £ 
(1-7,1/2), j=l,.-,r 



AA z j)-= SU P W v Z a Zj){u)\ < 

0<u<v<t 



Thus, the integral with respect to Z ; can be defined by relation Q, and it admits the 
following estimate: 



[ f(s)dZj(s) 

J a 



<C a A b , a (Zj) 



a \(s-a) a J a (s-l) 



\a+l 



-dz ds. (4) 



1.2 Assumptions 

Throughout the paper, C will denote a generic constant which may change from line 
to line. To emphasize the dependence of the constants on some parameters, we will 
put hem into subscripts. By |-| we denote both absolute value and the Euclidian norm 
in by (•,•) - scalar product in K", for any n > 1. Also, throughout the paper 
a e (1 - 7, 1/2) will be fixed. 

We will assume that coefficients of Q satisfy the following hypotheses. 

(Ml) for all t e [0,T},xeR d ,k=l,...,m, j = l,...,r 

\a(t,x)\ + \h(t,x)\ + \cj(t,x)\ < C(l + |jc|); 

(M2) for all t e [0,7], x,y € M. d , i = 1,... ,d, k= 1,. . .,m, j = 1,.. . ,r 

\a(t,x) -a(t,y)\ + \b k (t,x) -b k (t,y)\ + \cj(t,x) -cj(t,y)\ 
+ 1 d Xi cj (t ,x) - d x .Cj (t ,y) | < C \x - y\ ; 

(M3) for some j8 € (1-7,1] and all t,s e [0,T], x e R d , i=l,...,d,k= l,...,m, 
j= 1, 

\a(t,x) — a(s,x) \ + \bk(t,x) — b k (s,x) \ + |c ; -(f ,x) — cj(s,x) | 
+ \d x .Cj(t,x) - d Xi Cj(s,x)\ < C\t-sf . 
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1.3 Unique solvability of mixed stochastic differential equation and limit theorem 

In order to formulate the results, we need to introduce the following norm for a vector- 
valued function / 

||/L,= sup f f S \f(s)-f(z)\(s-z)- l - a dz). 
sero.fl \ Jo J 



and a seminorm 



l/llo '' = o<:< p v<fl (v-»)i-« + l (z-ur« dz 



The first result is about existence and uniqueness of solution. 
Theorem 1 Equation |2]l has a solution such that 

||Z|| oor <oo a .s. (5) 

This solution is unique in the class of processes satisfying (BJ. 

Another result we need is a limit theorem. Consider a sequence of equations 

X"{t) =Xo + J (a{s 1 X n (s))ds + b(s,X"^))dW(s)+c{s,X"^))dZ"{s)yt € [0,7*], 

where {Z„,n > 1} is a sequence of almost surely y-Holder continuous processes. 

Theorem 2 Assume that \\Z — Z n \\ T — > in probability. ThenX"(t) —>X(t) in prob- 
ability uniformly in t. 

The proofs of Theorems [T] and [2] are exactly the same as those in one-dimensional 
case given in l20l . 



2 Main results 

2. 1 Stochastic viability of solution to mixed stochastic differential equation 

We remind the notion of stochastic viability: for a non-empty set D C W 1 , process 
X = {X(t),t>0} is called viable in D if P(X(t) eD,t>0) = 1 whenX(O) efl. 

Assume that D is smooth, i.e. there exists a function <p : W 1 — > M, (p eC 2 ' 
such that its gradient d x (p(x) ^ when <p(x) =0 and 

D = {x : <p(x) >0}. 

Denote dD = {x : (p(x) = 0}. 
Theorem 3 Under assumptions 
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(VM1) foranyt>0,xeD 

j m d 

((p'(x),a(t,x)) + - 52 E b ki (t,x)b k i{t,x)(pl l l (x) > 0; 

1 k=\H=\ 

(VM2) for any t > 0, x G D, k = 1 , . . . , m, j = 1, . . . , r 

((p'{x),b k (t,x)) = ((p'(x),Cj(t,x)) =0. 

Then X is viable in D. 

Proof Since H-X^y < °° a.s., the process X is a.s. continuous, so it is enough to 
prove that for all t G [0, T) P{X(t) GD) = 1. 

For x G R , « > 1, we denote k„(x) = jn(\x\ An), 

Z n (t)=n [ k„(Z s )ds 

J(t-l/n)VQ 

mdk„(Z) = {k n (Z t ),te[0,T}}. 

It was proved in [20, Lemma 2. 1] that 

\\Z"-k„(Z)\\ 0T <CK y {k n {Z))n l -y- a , (6) 

where Ky(g) = sup 0<1<r<r \g(t) — g(s)\ j(t — s)? is the Holder constant of g. Note that 

\\Z"-Z\\ 0T <\\Z"-k„(Z)\\ QT + \\Z-k„(Z)\\ QT . 

Since Z is a.s. continuous, it is bounded, and hence ||Z — k„(Z)\\ T —> a.s., n — s- °°. 
This relation together with (|6) and the inequality Ky{k n {Z)) < Ky(Z) < °° gives 

\\Z" -Z\\ QT ^0,n^oo 
a.s. Defining X" as a solution to 

X"{t) =X Q + J' (a(s,X n (s))ds + b(s,X''(s))dW(s)+c(s,X"(s))dZ n {s)^,t G [0,T], 

(7) 

and applying Theorem [2] we have that X n (t) — > X(t) in probability uniformly in 

re [o,t], n^-oo. 

Equation |7]) can be transformed to the Ito stochastic differential equation 

X"(t) = Xo + J (a n {sX'(s))ds + b(s 1 X"(s))dW(s) 1 

where 

r 

a n {t,x) = a(t,x) + 52 c j(t^)Zj(t) 
= a(r , jc) + 52 Cj (t , x)n (k nj (Z(f )) - A»y (Z( (f - 1 /n) V 0) )) , 
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and k n j(x) is the j'th coordinate of k n (x). 

Since k n (Z) is bounded, the coefficients of this equation satisfy assumptions 
(H1)-(H4) in Appendix. It follows from (VM1) and (VM2) that the conditions of 
Theorem [7] are satisfied. Hence, for all t G [0,T], we have P(X"(t) 6 D) = 1. Since 
X"(t) -» X(t) as n — > °° in probability and D is closed, we get P{X(t) ED) = 1, as 
required. 



2.2 Positivity and comparison theorem 

As the first application of Theorem[3] we consider conditions supplying positivity of 
solution to equation Q. This question is of particular interest in financial modeling, 
where the prices of most assets cannot become negative. 

Theorem 4 Assume that for some d' < d and any i = 1 , . . . , d' 

(PI) Xoi > 0. 

(P2) for any x G M. d such that xi = 0, Xi > 0, / = 1, . . . , d', it holds a, (f , *) > and for 
any k — l,...,m, j — 1, ...,r b^t^x) — cp(t,x) = 0. 

Then P(X f (t) > Q,i = 1, . . . ,d',t e [0,T]) = 1. 
Proof Consider the equation 

X(t)=Xo+ / (d(s,X(s))ds+£b k (s,X(s))dW k {s) + Y,Cj(s,X(s))dZj{s)), 

J ° V /c=l j=l ' 

t e [0,7], 



(8) 



where 



a(t,x) = a(t, \xi\ , . . . , \x d >\ ,x d ' +l , . . . ,x d ) , 
b(t,x) = bit , |*i | , ... , \x d >\ ,x d ' +l , . . . ,x d ) , 
c(t,x) =c(t,\xi\,..., \x d i | ,x d > +1 ,x d ) . 

It follows from Theorem|3]that the unique solution X to equation |8]l is viable in each 
D t = {xi>0},i = l,...,d'. 

Consequently, under (PI), we have P(Z,(f) > 0,f G [0,T]) = 1 for i = l,...,d', 
or, equivalently, P(X,-(f) > 0, i = 1,„ . ,d',t G [0,T]) = l.Hence,X(f) solves and, 
by uniqueness, we have P(X(t) =X(t),t G [0, !T]) = 1, concluding the proof. 

Now we turn to a comparison theorem for solutions of stochastic differential 
equations. 

Let I G {1 , . . . , d} be a fixed coordinate number and X q , q = 1 , 2, be solutions to 
mixed equations 

X«(t)=X* + / (a«(s,X<i(s))ds+ £ b k (s,Xf(s))dW k (s) + £ cy^X^))^)), 



(9) 



8 



Alexander Melnikov et al. 



The coefficients of both equations are assumed to satisfy hypotheses from Subsection 
1.2. 

We note that the coefficients b and c are the same for both equations and they 
depend on the Zth coordinate only. This assumption seems restrictive, and we discuss 
it in Remark Q] 

The following result can be deduced from Theorem [3] exactly the same way as 
Theorem [8] is deduced from Theorem [7] 

Theorem 5 Assume that 

(CM1) X*,<X*. 

(CM2) for any x 1 ,x 2 £ W 1 such that x) — x 2 it holds a 1 ^^ 1 ) < a 2 (t,x 2 ). 
Then P{X}(t)<X 2 {t),t £ [0,T]) = 1. 

3 Applications 

To apply a pathwise comparison theorem to the mixed financial market model, con- 
sider the following one-dimensional pure and mixed analogs of the Cox-Ingersoll- 
Ross market price model, i.e. the stochastic differential equations of the form 

dX(t) = aX(t)dt + oX(t) x dB H (t),t >Q,X(Q) =X >0 (10) 

and 

dX(t) =aX {t)dt + aX(t) x (dW{t)+dB H {t)), t >0,X(0) =X Q >0 (11) 

with 1/2 < A < 1. The diffusion coefficient b(s,x) — c(s,x) — Ox x in both equations 
is not differentiable at zero and therefore does not satisfy (M2). Therefore, we can 
not apply Theorem 1.1 as well as any other known results concerning the existence 
and uniqueness of the solution of mixed stochastic differential equations involving 
fractional Brownian motion, to (jTOJ and ( fTT| . 

Suppose for the moment that equation (jTOJ has a solution X. Consider the integral 
JqX(s) dB H (s). According to the results concerning the pathwise Holder properties 
of the integrals with respect to ffim (for example, see [9]), this integral is Holder 
continuous in t up to order H. It follows that the trajectories of X have the same 
property. So, the process X x has Holder continuous trajectories of order up to HX. 
In turn, the integral [QX(s)^dB H (s) exists as a path-wise Riemann-Stieltjes integral 
if the sum of the Holder exponents of the integrand and the integrator exceeds 1 (see, 
e.g., [26|), i.e. if HX +H > 1. This heuristical argument explains why we consider 
equation ( fT0) i only for the values of index H G (yiXi !)■ Similarly, the solution of 
equation ( fTT| has Holder continuous trajectories up to order 1 /2, and the process 
has Holder continuous trajectories up to order ^, and hence, it must be 4 + H > 1, 
or H e (1 - |, 1). In particular, if X — 1 /2, then H G (2/3, 1) for equation {TO), and 
H G (3/4, 1) for equation (jll}. 

Consider ( fTT| on some interval [0, T] assuming additionally that the processes 
W and B H are independent. According to 0, there exists a Wiener process W with 
respect to the filtration generated by the sum W +B H , such that 

W{t)+B H {t) =W(t)- f f r(s,u)dW(u)ds, (12) 
Jo Jo 
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where the square integrable kernel {r(f,s),0 < s < t < T} is the unique solution of 
the equation 

r(t,s)+ [ r(t,x)r(s,x)dx = H(2H-l)(t-s) 2H - 2 1 < s < t < T. (13) 
Jo 

In this case, equation ( fTTj ) can be reduced to the stochastic differential equation 

dX(t)= (aX(t)-oX(t) x J r(t,u)dW(u)jdt + oX(t) x dW(t), ^ 

t >0 } X{0) =X >0. 

Equation ( fl4] > does not involve a fractional Brownian motion but it has a random 
non-Lipschitz and non-Markov drift coefficient. 

Theorem 6 (I) Equation (|T0]> has a unique solution ifH € (yxji !)• 
(II) Equation \\\) has a unique solution ifH £ (1 — 4,1). 

(Ill) Denote Vo = inf{f > : X{t ) = 0}, where X is a solution to either ( |T0] > or ( fTTj i. 
77ie« X(f ) = a.s. for all t > Vo- 

Proof The proof of this theorem is given in Appendix C. 

Consider the one-dimensional mixed Cox-Ingersoll-Ross interest rate model de- 
scribed by equation ( [TTj ). We suppose that the discounting factor equals B(t ) = e at , t > 
0. It is impossible to solve explicitly equation ([TTJ and the distribution of its solution 
has a complicated form. We note that even for the classical Cox-Ingersoll-Ross inter- 
est rate model, involving only Wiener process, it has a noncentral chi-square distri- 
bution according to [6| (see also a discussion in O). To avoid such difficulties, we 
apply Theorems|6]and|5]to get an upper bound for option prices on interest rate in the 
model described by equation ([TTJ. To this end, introduce Y(t) =X(t) 1 ,t € [0, Vo). 
According to the Ito formula, Y satisfies the following stochastic differential equation 
on [0,V ): 

dY(t) = (fl(l - l)Y(t) - ^^(t) X) ) dt + °( 1 -V{ d W(t)+dB H (t)). (15) 

Consider an auxiliary process Z that satisfies the following stochastic differential 
equation (Vasicek model) 

dZ(t) = a(l -X)Z{t)dt + (7(1 -X)(dW(t) + dB H '(f)) . (16) 

with the initial condition Z(0) = Xq , Equation ( fT6"J has a unique solution on R+, 
and it is a Gaussian process of the form 

Z(f )=X*- X e ail - l), + a(l- A) J\ a{l - l){t - s] (dW{s) +dB H (s)). 

Applying Theorem[5]to Y and Z on the interval [0, Vo), we get Y(t ) < Z(t) a.s. for all 
t G [0, V ). Since Y{t ) = for f > v , this implies Y(t) < Z{t) + for all f > 0. Therefore, 
X(t) < Z(t) V^ 1 ^ a.s. for f > 0. Considering a European option with a nondecreasing 
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payoff function /( •) and maturity T > 0, we can use the latter inequality to derive the 
following upper price bound: 

E [e-" T f{X T )} < E 
where the distribution of Z(T) is Gaussian with mean 

X i-X e a(\-X)T and 

variance 

g 2 (l-A) 2 ( \ - ~ +H(2H-l) / e« l -W'+^\t-s\ 2H - 2 dtds). 

\ 2a(l — A) Jo Jo J 

For a = 0.1, Xq = 1,A = 0.5, T = 10, f(x) — (x — K) + we have calculated the 
upper bound for different a and K. We have also computed numerically the prices 
by using the Euler method with 4096 points to simulate 20000 paths of the solution. 
Results are following. 



a 


0.1 


0.5 


1 


K 


0.5 


1 


2 


0.5 


1 


2 


0.5 


1 


2 


Price 


0.8818 


0.7015 


0.4032 


2.17 


2.04 


1.891 


4.448 


4.42 


4.312 


Upper bound 


0.8953 


0.7198 


0.4192 


2.55 


2.434 


2.223 


6.552 


6.448 


6.25 



As we see, the difference between the upper bound and the price is minor for 
small values of a and becomes substantial when a increases. There are two reasons 
for that: first, the term we have omitted in ( fT5j ) to obtain the upper bound grows with 
a, second, for large a the solution hits zero quickly and stays there, while the upper 
bound becomes larger. 

It is of great interest whether or not solutions of equations ( fTT| and ( fTO) stay pos- 
itive. For equation (JTTJ with H > 3/4 a complete answer can be given. Namely, in 
this case the sum W + B H under equivalent measure trasformation becomes a stan- 
dard Wiener process, so equation ( [TT| transforms to an Ito equation. Its solution is a 
diffusion process, so it can be checked by using the results of lfT3l Section VI. 3] that 
the solution almost surely vanishes within a finite time. 

For equation ( fTO) , the precise answer is not known. The numerical experiments 
below suggest that for a < the solution vanishes almost surely, and for a > it 
stays positive at least with positive probability. (The Euler method we have used for 
simulations cannot give and answer whether the probability is 1 or less.) 

We have considered equations 

dX, = aX,dt + VX~,dB^,X = 1 , t > 0, 

with H = 0.8 and a = ±1. For a = 0.1, the following figures contain histogram of the 
time when trajectories hit on the segments [0,50] and [0,500]. 



e- aT f(Z{T)f l - X) ) 
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□ 5 10 15 20 25 30 35 40 45 50 50 100 150 200 250 300 350 400 450 500 



We see that almost the same number of paths (around 400) stays positive at t = 50 
and at t = 500. 

For a = —0.1, the following figures contain the graph of the 1000 paths and his- 
togram of the time when trajectories hit 0. 




50 100 150 200 250 300 350 400 450 500 50 100 150 200 250 300 350 400 450 500 



We see that only a signle path stays positive at t = 500. 



A Stochastic viability for ordinary stochastic differential equations with 
random coefficients 

Consider a stochastic differential equation 

Xit) =X Q + [' (a(s,X(s),(o)ds+ V b k (s,X(s),w)dWk(sj), t e R + (17) 
Jo \ fcf] ' > 

where the coefficients a and foj-, k = 1. . . . ,m, are measurable functions from R + X X £2 to R d ; W^, 
k = 1, . . . ,m are independent F- Wiener processes, Xq is J?o-measurable random vector in W 1 . We assume 
for simplicity that Xq is bounded. 

The following assumptions guarantee that there exists a unique solution to \17\ in the class of adapted 
square integrable processes (see |22|). 

(HI) For each t 6 1R + the functions a(t, ■) and £>t(r, ■), k = 1. . . . ,m, are .^(W 1 ) <g> ^-measurable. 
(H2) Functions a(-,w) and ^(-, to), k = 1. . . . ,m, are almost surely jointly continuous on K + x W 1 . 
(H3) For any t G K + and any x e R 

|a(t,0)| + £|Mf.0)l<C. 

k=l 
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(H4) FotmyteM. + ,x,yeW 

m 

\a{t,x)-a{t,y)\ + Y,\h{t,x)-b k (t,y)\<C\x-y\. 

k=l 

Let <p : R d -> M, <p G C 2 (R d ) be a function such that its gradient d x q>(x) / when <p(x) = 0. Let the 

set 

D = {x: (p(x) > 0} 
be non-empty, and dD = {x : (p(x) = 0} be its boundary. 

Theorem 7 Assume ifta/ Xq E D a.s. and the following conditions hold almost surely: 
(VI) for any t E R+, x E D 

| in d 

a(t.x) := {(?'(x).a(t,x)) + - £ £ h ki (t,x)b H {t,x)<f^{x) > 0; 

Z £=1 ,-,(=1 

(V2) /or any r e IR+, x€D,k=l,...,d 

P k (t,x):=(<p'(x),b k (t,x))=0- 

Then P(X(t) E Dfor allteR + ) = 1. 

Proo/ Since D is closed and X(f) is a.s. continuous, it is enough to prove that P(X(t) E D) = 1 for all 
t>0. 

Step 1. First we prove that if a(t,x) > a.s. for all t E R+, x E D and if Xo G 3D, then there exists 
a stopping time 9 > a.s. such that X(t) G Z) for all t G [0, 6]. To this end, fix some R > \Xq\ and define 
stopping times 

T = inf {s G R + : ce(.v,X(.s)) < 0} , z R = min{.v G R + : |X(s)| > R} . 

As usual, we suppose that a stopping time equals ~ if the corresponding set is empty. It is clear that z > 
and z R > a.s. 

For any u > put Q u = u A z A z R and apply the Ito formula to the process <p(X(-)): 

<p(X(6 u )) = J o (a(.v,X(.v))rf.v+£j3,(.v,X(.v))dW i (.v)j. 

Since X is bounded on [0, 9 U ] , the expectation of the stochastic integral equals zero. Hence 

El<p(X(t))] = E\j\( S ,X(s))ds . 

For a non-negative function y G C(R) such that J R y(x)dx = 1 and y(x) = 0,x^ [0, 1], define 

\)/„(x) = n /' ' / y(nz)dzdy. 
Jo Jo 

Obviously, ift(i) "t" as « — ^ <=° and \y' n (x)\ < 1, n > 1. 
Applying the Ito formula to tfr„(<p(X (■))), we get 

r% m 
2 ./o t f i 
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Similarly, 



E[ ¥n (<p(X(9 u )))] = E 
1 



+ - 



[ y n (.<p(X( S )))a(s,X(s))d S 
Jo 

£E / ^'(<p(Z(,)))i3,(.v,X(.v)) : 

Jo 



(18) 



Recall that a (s, X s ) > for 5 < 9„ , and (x) < 1 , so the first term in the right-hand side of (18} does not 

exceed E f®" a(s,X(s))ds\ . We will prove now that the second term vanishes. 

For x 6 Mr , let z(x) be the closest to x point such that (p(z(x)) = (any of the points if there are more 
than one). As the function <p is continuously differentiable and its derivative is non-zero on dD, there exists 
e > such that 

\x - z(x) I < C R \<p(x) - <p(z(x))\ = C R \<p(x) I 
whenever \x\ < R, \ (p(x)\ < e. Further, for \x\ < R 

\p k (s,x)\ = \Pk(s>x)-P k (s,z(x))\ 
< I {<p'(x) - <p'(z{x)),b k (s,x)) I + I {<p'{z(x))Ms,x) - b k {s,z(x))) 
< Cr(\x-z.(x)\ + \x-z.(x)\) < C R \x—z(x)\ . 



But for n > 1 /e it holds that 



whenever \(p(X(s))\ > e, whence 



irf((p(X(s)) = nyr(n<p(X(s)))=0 



f 9 " lrf(<p(X(s)))P k ( S ,X( S )) z } ds < C R n f" ¥ (n<p(X(s)MX(s)) 2 ds 
Jo Jo 



. CrU 2 

< sup a: mx) -> 0. n— > <*>. 



Letting n — > °°, we get from j 1 8| that 

E[\(p(X(e u ))\] < liminfE[v/„(9(Z(e„)))] < E 



/ a(s,X(s))ds 
Jo 



= E[<p(x(e u ))], 



and hence, <p(X(6 u )) > a.s. Since 0„ = u A T A T« and « is arbitrary, we get the desired claim with 
8 = TM/j>0 a.s. 

Step 2. We prove the statement of the theorem under assumption that a(t,x) > a.s. for all t > 0, 

x e D. 

Define To = inf{.s > : X s ^ D} and assume on the contrary that P(To < °°) > 0. Then for some 
R > it holds P(t d < °°,|X(t )| < R) > 0. Consider a new stochastic basis (i2,^",F,P), where £2 = 
{t < oo, X(t d ) < R}; ^ = F = {^"r+T D |fl,f > 0}, where ^" r+Tz , is the cr-algebra generated by the 
stopping time t + T D ; P(A) = P(A)/P(.fi), A 6 Define also X(r) = X(t + T D ), a(f,;t) = a(r + T ,jt), 
b k (t,x) = b k (t + T D ,x),k=l,...,m, W(t) = W(t + T ) - Evidently, W is a F-Wiener process, and 

the newly defined coefficients satisfy (H1)-(H4), (VI), and (V2), moreover, a(t, x) > P-a.s. for all t > 0, 
jt 6 D. It is easy to see that X solves 

X{t)=X(0)+ [a(s,X(s),w)ds+Tb k (s,X(s),m)dWk(s)), t>0. 
Jo V t t5 ' 

According to Step 1, there exists 6 > P-a.s. such that E D, t e [0,9] P-a.s. Thus, X(f) 6 D, 
f 6 [To, T/j + 0] for P-a.a. CO 6 12, which contradicts the definition of 'Co- 
Step 3. Now we prove the statement in its original form. Let {a"(t,x),n > 1} be sequence of coeffi- 
cients such that (H1)-(H4) holds, for any T > 



sup sup \a"(t,x) — a(t,x)\ 
te[o,T\x, y m 



•0, n- 
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and almost surely for any t >0, x E D the assumption of Step 1 is satisfied, i.e. 

I m d 

a"(t,x) := (<p'(x),a«(t,x)) + - £ £ b a {t,x)b Kt (f,x)tf } (x) > 0. 
*=lij=l 

One can take, for example, a n (t,x) = a{t,x) + ra -1 (p'(x)G(x) with a positive smooth function G: R rf — > R 
which does not vanish on dD and decays on infinity sufficiently rapidly so that (p'(x)G(x) is bounded 
together with its derivative. 
Let X" be the solution of 

ft i "' \ 

X»(t)=X Q + U{ S ,X"(s))ds+Yh(s,X n {s))dW k {s)), t>0. (19) 
Jo V £ I 

Then it is well known that X" (f) X(t ) in probability locally uniformly in t. By Steps 1 and 2, X" (t) e D 
a.s. for all f > 0. Since D is closed, we get X(t) 6 D a.s. for all f > 0. 

We note that one can extend the results of Theorem[7]to the case where the constants in (H3), (H4) are 
random and (H4) holds locally in x,y. 



B Comparison theorem for equations with random coefficients 

We formulate the result below as multidimensional, while it is basically about a pathwise comparison in 
one-dimensional case. 

LetX 9 , q = 1,2, be solutions of stochastic differential equations 

ft i "' \ 

X«(f) = Xg + J q [tfi(i,Xt(s))di + £ b k ( S ,X?(s))dW k (s)) , r > 0, 

where the coefficients a*, q = l,2, wdb k , k = l,...,m, satisfy (H1)-(H4); X$ = (X^, . . . ,X* d ), q = 1,2, 
are bounded .^-measurable random vectors, and I £ {1, ...,d} is a fixed coordinate number. 

Theorem 8 Assume f/iar 
(CI) X^<X 2 a.s., 

(C2) /or any x',jc 2 e R <i such thatxj =xj it holds a [ (t,x [ ) < a 2 (t,x 2 ) a.s. 
Then P(X/(t) < Xj(t),t > 0) = 1. 

Proof Consider the process X(t) = (X 1 (t),X 2 {t)) € R M and set <p(x) = xj -xj. Then, in the notation 
of TheoremF] D = {X, 1 < X 2 }, dD - {Xj = X 2 } and a(t,x) = a 2 (t,x) - a 1 (t,x). By the assumption, 
Xq = (X l ,Xq) 6 D and (x(t,x) > a.s. for x 6 5D. Thus, we get the desired statement from Theorem]?] 

Remark I In 1 18], the following result is given (we reformulate it slightly according to our notation). 

Theorem 9 (| 18 1) Let X q ' , q= 1,2, be solutions of stochastic differential equations 

ft , \ 
X"(t) = X*+ / U<(sXi(s))ds+ £ bl(s,X"(s))dW k (s)) , t > 0, 
JO v k=[ ' 

where a q : R + x W 1 — > R rf and b q : R + x W 1 — > R'' satisfy the linear growth and the Lipschitz continuity 
assumptions, X ' 6 M are non-random. Let I C { 1 , . . . ,d} be a non-empty set. Then the following assertions 
are equivalent: 

(a) for any X* = (Xg,, . . . ,X^), q = 1, 2, swcft dial X^ < X^/or all jel it holds 

P(X}(t)<xj(t),jeft>0) = l; 
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(b) for any p E I and x' = (x\, . . . ,x' d ), i = 1, 2, such that < x~j, j G /, x p = x 2 it holds 

a),(t,x l )<a 2 p (t,x 2 ), b [ kp (t,x i )=b 2 kp (t,x 2 ),k=i,...,m. 

First we note that it follows from (b) that b' kp , i= 1,2 are equal and independent of all coordinates except 
those from /. 

Now we argue why comparing more than one coordinate makes only little sense. For example, take 
/ = { 1 , 2} and observe that b' kp are independent of x. Indeed, from (b), for any p G / 

>>lp(t,xi,X2,...) = bl p (t,x l ,x 2 V4i---) =b\p(f,x\,xi V4i--.) 
= b\ p {x y Vx' 2 ,x 2 Vx' 2 ,...) = bl p ( Xl Vx' 2 ,x 2 Vx' 2 ,...), 

and similarly 

b l kp {t,x u x 2 ,...)=b l kp {x l Vx' 2 ,x 2 Vx 2 ,...). 

The dots here can be anything, so we derive b\ (t,x) = b\ „{t,x) = b kp (t). Thus, for /)£/, defining Y'(t) = 
X' p (t) — Bpit), where B p (t) = Y!'k=\ Jo bkp{s)dWk{s), we can write the equation for Y' as 

dY'(t) = a'(t,Y'(t)-B„(t))=:f(t,Y'(t))dt. 

Though /' is random, this is an ordinary differential equation. We can solve assuming that all other co- 
ordinates and a are given. The conditions we have guarantee that f l (t,y) < f 2 (t,y) disregarding what 
happens to other coordinates or CO, so from comparison theorem for ordinary differential equations we 
conclude that Y { (t) < Y 2 (t) for all t > whenever F'(0) =X^ < Y 2 (0) = X% p , whence the same is true 
forX;(f). 

Now for / containing a single element we stillget that b(t,x) does not depend on other coordinates. For 
a, we get exactly condition (C2) from Theorem Isl Again, assuming that all other coordinates of X 1 and 
X 2 are given, we get one-dimensional equations with adapted random coefficients. But (C2) guarantees 
that this random coefficients are properly ordered whatever happens with other coordinates, so this is 
essentially a one-dimensional result in a very rigorous sense. It follows that there is nothing new in this 
result compared to [12] (except randomness of coefficients, but it can be checked that the proof in 1121 
works for random coefficients). 



C Proof of Theorem© 

To avoid technical difficulties, we consider equation (10) with X = j . The proof is similar for remaining 
cases. 

Cases (I), (II). For any e G (0,Xo) we introduce a smooth coefficient of the form 

d e {x) = sfr\ x >e + (5/2£- 3 / 2 .v 2 - 3/2e- 5 ' 2 x 3 )\o< x < e - 

The function d e G C 1 (M), and its derivative 

d' E (x) = l/2x- l ' 2 h> £ + (5e- 3 ' 2 x-9/2e-^ 2 x 2 )l < x < £ 

satisfies the Lipschitz condition. Therefore, the following equation with the smooth diffusion 

dX £ (t) = aX £ (t)dt + d e (X e {t))dB H (t),t > 0,X £ (0)=X > (20) 

has a unique solution. Let us introduce the Markov moment 

T £ = irjf{/ > : X e (t) < e} > 0. 

Obviously, X e (t E ) = e on the set {t £ < °°} and X e (t) > e for t > on the set {t e = <»}. It means that on 
the interval [0, T E ) there exists the solution of equation \ W) because on this interval coefficients of JTo} and 
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420| coincide. Moreover, the solution of (10) is unique on [0, T E ) because any solution of on [0, T e ) 
can be extended to the solution of j20} on M + . In turn, it means that solution X of equation \10) exists and 
is unique on the interval [0, sup e>() T e ). We note also that T E] < T E2 for £2 < £1. In order to prove that X 
is continuous at the point sup e>0 T e on the set {sup E>0 T e < °°} (in this case X(sup E>0 T e ) = on the set 
{sup E>0 T e < 00}), we estimate the expectation of the solution. 

Let #> 2/3, e>0,f >0,R> 1, L> and J3 e (1 -H, 1/3) be fixed. For 5 e (0, 1/3 -j8), consider 
the process 



Jo it 



X(t)-X(u)\ 



(iu. 



It is well defined on any [0, T e ] because the numerator does not exceed C((Q)(t — u) H K for any < K < H 
and we can choose K in such a way that H — K > 2/3 + S. Define 

V L = wf{s >0:X(t) + Y(t) >i}AT f , 
v R = wf{s > : A sfi (B H ) > R} 

and 9 = 6r l = Vl A Vr. According to j4}, we have 



supZ(zA6) 

<X Q + \a\ / E[X(sA9)}ds + aE sup / ^/X[s)dB H (s) 
Jo l z < t Jo 

< Xo + \a\ J' E [X(s A 9)] ds + CpRo ( J' ^ E ^ Ad ^ ds 



The integral 1\ = J ' 



J5 



v / E[|X(iAe)-X(uAe)|] 



rfi admits an upper bound 



duds 



h<C f E[X(sA9)]ds + O l - 2IS . 

Jo 

To estimate the integral I 2 = Jg Jq V E \i x ^ B J-^ uAe ^ duds, note that 



v / E[\X(sA9)-X(uA9)\] 



(s-u) [ +P 
1 fE[\X(sA9)-X(uA9)\] 



+ (s-u)- [+s ). 



whence 

Therefore, denoting 
/iW = E 



h <C + 



{s- u y+ 2 i>+ 5 

E[\X(sA9)-X(uA9)\] 



Jo Jo (s-uy+ 2 P+ s 



duds. 



supX(zA0) 



f 2 {t) = E[Y{tA9)}= f 
Jo 



E[\X(tA9)-X(uA9)\] 



{s-uy+w+s 



du. 



we obtain from [21) and the above estimates that 

f\(t)<C + CR [' fi(s)ds + CR f ' f 2 (s)ds 
Jo Jo 



(21) 



(22) 



(23) 
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Furthermore, the function f% admits the following upper bound: 



lulu {s-rY+P J 

f ' M S ){t-s)- 2 P- s ds+CR f /M/) S -^(t-s)-^- S ds 
10 Jo 



(24) 



+CR f r^MM^Ewm dr ^ s) -^ ds . 

Jo Jo • " 



Consider the second and the third terms in the right-hand side of {24j separately. Recall that 5 < 1/2 — /3 
and rewrite the second term as follows: 

Jo 

Jo 

= ['Ms)( t -s) i - 6 p- f > s d S + ['s-W(t-sr H 

Jo Jo 

= [' h(s){t-sf-^- 6S ds + Ct 2S . 
Jo 

From the choice of /3 we have 1 — 6/3 — 65 > — 1 . To estimate the third term, we apply J22) and get 



(25) 



ds 



Jo Jo (s-r) l+ P 

< f f 2 M(t- S )- 2 P- s d S + Ct [ - 2 P. 
Jo 

Finally, we get from j2lj-j26j that 

/iW+/2(f)<CT?(l + £(M S )ds+f 2 (s))(t-syds), 



(26) 



(27) 



where 7= (1 — 6/3 — 65) A (—2/3 — 5) > — 1. It follows from \21\ and the generalized Gronwall inequality 
that 

ft W +/2W < C#exp{Cft}. (28) 
Furthermore, the upper bound in J28l does not depend on L. Therefore, by the Fatou lemma, 



supZ(zA Vr AT £ 



and by the Lebesgue monotone convergence theorem, 



< CRexp{CRt}, 



supZ(zAv«A(supT £ )) 



<CRexp{CRt}. 



It means that the both components, fyXfyds and J f y / X(s)dB H (s), are almost surely bounded and con- 
tinuous on [0, V/; A (sup E>0 T e )] if sup E>0 T E < 00 • Since Vr — > °° a.s., we claim that X is continuous on 
[0, sup e>0 T E ] if sup e>0 T e < °°. But X(t e ) — > as e — > if sup E>0 T E < °°, therefore X(sup e>0 T f ) = on 
the set {sup e>0 T e < °°}. 

Case (III). It is easy to see that Vo = sup E>0 T e . Assume that there exists the solution of equation ACT 
with X = i such that for some p > probability of the event A = { Vp = inf{f > Vo : X(t = p} < °°} is 
nonzero. Then we can consider the new probability space Q = A and repeat the arguments above to prove 
the uniqueness of the solution of equation \ IQ) on some interval [v p , v). However, equation (10) has on 
this interval zero solution, and we get a contradiction. 
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